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ABSTRACT 

We present results from a study of short term variability in 19 archival observations by XMM- 
Newton of 16 Ultraluminous X-ray Sources (ULXs). Eight observations (six sources) showed 
intrinsic variability with power spectra in the form of either a power law or broken power 
law-like continuum and in some cases quasi-periodic oscillations (QPOs). The remaining ob- 
servations were used to place upper limits on the strength of possible variability hidden within. 
Seven observations (seven sources) yielded upper limits comparable to, or higher than, the val- 
ues measured from those observations with detectable variations. These represented the seven 
faintest sources all with f x < 3 x 10~ 12 erg cirT 2 s _1 . In contrast there are four observations 
(three sources) that gave upper limits significantly lower than both the values measured from 
the ULX observations with detectable variations, and the values expected by comparison with 
luminous Galactic black hole X-ray binaries (BHBs) and Active Galactic Nuclei (AGN) in 
the observed frequency bandpass (10 3 - 1 Hz). This is the case irrespective of whether one 
assumes characteristic frequencies appropriate for a stellar mass (10 M Q ) or an intermediate 
mass (1000 M Q ) black hole, and means that in some ULXs the variability is significantly sup- 
pressed compared to bright BHBs and AGN. We discuss ways to account for this unusual 
suppression in terms of both observational and intrinsic effects and whether these solutions 
are supported by our results. 
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, 1 INTRODUCTION 

Ultraluminous X-ray sources (ULXs), point-like sources with lu- 
minosities greater tha n 10 39 erg s -1 , were first identified by the 
Einstein Observatory (Fabbiano 1989). Initial studies of short term 
variability in dicated that they were likely to be accreting com- 
pact objects I Ok ada et al. 1998). Both spectral analysis and long 
and short term variability studies suggested similarities between 
ULXs and Galactic black hole X-ray binaries (BHBs) (compris- 
ing a stellar mass black hole accreting matter from a p artner star 
dMiller & Coibertl |2004 lFabbiandl200d lRobertsll2007h ). The X- 
ray spectra of ULXs can be fitted by models commonly used 
for BHBs jColbert & Mushotzkvl 1 19991) . and display changes in 
their energy spectra si milar to the state changes observed in BHBs 
dLa Parola et "al] |200lh . The high luminosities are greater than that 
predicted by applying the Eddington limit allowed by spherical ac- 
cretion onto a typical stellar mass black hole (< 20 M Q ). ULXs are, 
however, mostly observed outside the nucleus of their host galaxies 
and thus are unlikely to be super massive black holes (> 10 5 M Q ), 
typically observed in the centre of galaxies as active galactic nu- 
clei (AGN). Their high observed luminosities and the low appar- 
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ent disc temperatures measured from the energy spectra have led 
to suggestions that these sources represent an until-recently un- 
known class of black holes with masses between those of stellar- 
mass black holes and AGN around 100 - 10000 M : inte rmediate 
mass black holes (IMBHs) fcolbert & Mushotzkvl [l999h . In con- 
trast a number of theories have been put forward suggesting that 
the high luminosities may be caused by a stellar mass black hole 
emitting or accreting i n different ways. Emission processes such as 
anisot ropic emission JKing et al. 2001), relativistic beaming from 
jets jGeorganopo ulos et al. 2002), or accretion onto the BH over 
the expected limit (Beselman 2002) could explain some observed 
features. Constraining the masses of these objects would be an im- 
portant step forward in understanding them. 

A number of different relationships between BH mass, ac- 
cretion rates and features of power spectral densities (PSDs) 
common to both BHBs and AGN have already been identified 
( McH ardv et al.l2006l;lMcClintock & Remillard2006l;lyan der Klis 
20061; [Done & Gie rlihski 2005; Marko witz et alj2003l ; IUttlev et al. 
2002; Belloni & Hasinger 1990). If these key features are seen in 
the PSDs of ULXs then it may provide a way to estimate the black 
hole mass. BHBs are observed in a number of different "states" tra- 
ditionally defined by the dominance of components within the en- 
ergy spectra, which also show distinctly different PSDs. It is now 
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recognised however that the accret ion rate may in fact be a more ef - 
fective parameter to discern them dMcClintock & Remillardl2 006). 
Hardness-Intensity diagrams following BHBs through outbursts 
have revealed a form of hysteresis p ath that these sources follow 
as they move between different states dMivamoto et al.|[l995h . 

A typical outburst cycle begins in a 'low/hard' state in which 
the X-ray spectrum is dominated by a hard power law, with a weak 
contribution from a thermal component originating in the accre- 
tion disc. The X-ray PSD shows band-limited noise that may be 
approximate d by a doubly broken power law or multipl e broad 
Lorentzians (Nowak 200d lMcClintock & Reimilardll2006l) . In this 
state the PSD often displays quasi-periodic oscillations (QPOs). 
Radio emission from a steady jet is frequently observed. At the 
peak of outburst the source may reach a 'very high state' in which 
the X-ray spectrum shows a steep power law and a strong thermal 
component. The PSD typically features a complex, broad contin- 
uum and strong QPOs. The third main state is the 'high/soft' state 
characterised by a dominant soft, thermal spectrum. In this state the 
PSD appears as a broken or bending power law with an index of -2 
above the b reak and -1 bel ow it, which may extend to very low 
frequencies dReig et al]|2003l) . QPOs are not usually seen and the 
radio emission is greatly reduced. Towards the end of outburst the 
source returns to the 'low/hard' state and the luminosity decreases 
until the source reaches quiescence. 

Observations of ULXs, particularly with XMM-Newton, have 
allowed for X-ray spectral and timing analysis which appears to 
show some of the features observed in BHBs. The PSDs from 
XMM-Newton are most sensitive to frequencies around which we 
would expect to see a high frequency break (present in both the 
high/soft and low/hard states) if the mass of the black hole is suffi- 
ciently large (i.e. the IMBH model predicts a break at ~ 1(T 4 - 10~ 2 
Hz). 

Observations in all three states of BHBs, as well as all well- 
constrained AGN data, show a high-frequency break above which 
the PSD decays as a steep power law (index < -2). This break fre- 
quency may be a key diagnostic feature in establishing the proper- 
ties of the system. McHardv et al. (2006) showed that the break fre- 
quencies in AGN and BHBs display a correlation dependent on the 
mass of the black hole and its accretion rate T Break cc M^ni^J^ . 
BHBs in the low/hard state have a sec ond lower frequency break 
from a flat power law to an index of - l. lKording et al.l j2007h h ave 
altered and extended the analysis from iMcHardv et al.l feOOtj) to 
BHBs in the low/hard state through modelling the PSD as a se- 
ries of Lorentzians (apparent when the PSD y-axis is plotted as 
P x f) and then forming a relationship between accretion, BH mass 
and the Lorentzian relating to the high frequency break. QPOs 
have been detected in two ULXs, M82 X-l and NGC 5408 X- 
1, indicating th at the systems may have similar variability prop- 
erties to BHBs dStrohmaver & Mu shotzkv 2003; Stro hmaver et al.l 
1 20071) . By contrast, other studies of ULX variability have indicated 
that these sources often do not show strong short term variability 
dFeng & Kaarell2005l : lGoad et alj2006l ; lRobertsll2007l) . However a 
break in the PSDs of ULXs NGC 5408 X-l and M82 X-l, from a 
flat power law to -1 dependence, similar to that seen in t he low/hard 
state o f BH Bs, has been identified b e low the QPOs bylSoria et al.l 
d2004l) and iDewangan etail d2006l) ; iMucciarelli et alj d2006h k- 
spectively. Identifying the higher characteristic break frequency in 
the PSDs of ULXs could allow for a more accurate determination 
of their masses, and a deeper understanding of the nature of these 
particular sources. 

This paper describes a power spectrum analysis of the time 
series obtained from 19 XMM-Newton observations of ULXs. Sev- 



eral of these observations have been analysed and discussed sep- 
arately elsewhere (see section 2), but this paper presents the first 
uniform analysis and presentation of the power spectra from a 
reasonable-sized sample of good quality ULX observations (i.e. 
observations with good exposure time of sources with high X-ray 
fluxes). As noted above, some ULX observations show no evidence 
for rapid X-ray varia bility, but in almost all published studies (e.g. 
iFeng & Kaaretl }2005) it is not clear whether the lack of detected 
variations is due to insufficient data or an intrinsically low variabil- 
ity amplitude (by comparison with well studied BHBs and AGN, 
and also t he variable ULXs ). Indeed, the careful analysis of Ho 
II X-l bv lGoadetaf]d2006t) showed the ULX to be intrinsically 
under-variable. This raises some interesting questions, but it is not 
clear whether the lack of strong variability is common in ULXs: in 
all other previous studies where no variability was detected, there 
was no attempt to constrain the amplitude of intrinsic variation al- 
lowed. In the present study we analysed the power spectra of all 19 
observations, and whe re no variab i lity w as detected, we followed 
the approach used by I Goad et al] d2006T) to constrain the power 
spectrum. This allowed us to compare the variability properties of 
all the ULXs, whether variability could be detected or not, in a way 
not possible using the results of previously published studies. 

The rest of this paper is organised as follows. Section 2 de- 
scribes our sample selection criteria and search process. In Section 
3 we discuss the basic data reduction techniques used and methods 
of temporal analysis. Section 4 explains our results in some detail 
including investigations into the upper limits of variation hidden 
within the power spectra. Finally Section 5 gives a brief discussion 
and conclusion. 



2 SAMPLE SELECTION 

We have performed a search for bright ULXs with a flux previously 
observed to be greater than 0.5x 10~ 12 erg s _I . A sam ple was formed 
throug h use of sources previously iden tified in IStobbart et al " 
(2006 ), searching the ROS AT catalogues of iLiu & BregmarJ 1 200: 



andOu & Mirabel] d2005l) and by carrying out an archival study of 
papers studying ULXs using the ADS database (ADsQ). A search 
was then made for long observations (> 25 ks) of regions including 
these sources available by February 2008 in the XMM - Newton 
public data archive (XS^O) . The final sample consists of 19 obser- 
vations from 16 sources with the potential to provide useful timing 
data (see TableQJ. 

As mentioned above, several of these observations have been 
discussed elsewhere in the literature, and here we briefly summarise 
these reports. The observations of NGC 1313 X-l and X-2, NGC 
5204 X-l, NGC 4945 X-2, NGC 4861 ULX and NGC 253 PSX-2 
have no published discussion of their short-term variability. 

The observations of NGC 55 ULX w ere known to display 
interesting variability ISto bbart et alj[2004h . but this has not pre- 
viously been analyse d in detail. No variability was detected by 
IFeng & Kaaretl d2005l) in the observations of NGC 4395 X- l NGC 
3628 X-l, M83 ULX and NGC 2403 X-l. iBarnard et all d2007h 
detected no variability in NG C 4559 X-l (revising the previous 
report by [Cropper et alj|2004h . Only the remaining five observa- 
tions have detailed power spectral analyses published to date. The 
PSDs from the two observations of M82 X-l have been published 
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in lStrohmaver & Mushotzkvl j2003h and lDewangan et aljjioo^ re- 
spectively: these show a strong QPO and a break in the PSD. A 
similar QPO and spectral break has also been detected in the obser- 
vation of NGC 5408 X- 1 dStrohmaver et alj2007h . |Pewangan et al.l 
(2006) claimed a modestly significant QPO in the power spec- 
trum of the obser vation Ho IX X-l. Finally, as mentioned above, 
lGoadetalJ (2006) presented a detailed analysis of the lack of vari- 
ability present in the long observation of Ho II X- 1 . 



3 ANALYSIS 

Light-curves and energy spectra were extracted from the EPIC pn 
camera on XMM-Newton. In all observations the camera was op- 
erating in full-field mode excepting that of Ho IX X-l (obs Id. 
0200980101) which was taken in large window mode. The MOS 
cameras were not utilised due to their lower time resolution, mak- 
ing it harder to constrain the Poisson noise level during power spec- 
tral analysis. When the data were included they were found to add 
little to either the results or analysis. The data from each observa- 
tion was extracted and processed using the XMM-Newton SAS ver- 
sion 7.1.0. Data was filtered leaving only events with PATTERN < 

4 and FLAG==0. Events were extracted either within the optimised 
radius of the source, as evaluated by the SAS, or a radius of our 
own definition if other factors (such as the source sitting close to a 
chip edge or near another X-ray point source) made this unsuitable. 
A background light curve was taken from a rectangular region on 
the same chip and as close to the source as possible. 

Light-curves for the 0.3-10 keV band were extracted with a 
time resolution of 73.4 ms for all sources except for Ho IX X- 
1, the light curve of which had a time resolution of 96 ms ac- 
counting for the different camera operational mode. Most of the 
ULXs are near the aimpoint of the image and are clearly detected 
as separate point sources, with a few exceptions. Most sources 
were therefore extracted from a region centred on the co-ordinates 
given in Table Q] NGC 4395 X-l is close to the edge of the 
chip; M82 X-l and NGC 3628 X-l are close to the centre of 
their host galaxies and other bright X-ray sources, causing diffi- 
culties in exactly identifying the source whilst excluding all others 
given the resolution of XMM-Ne wton. For M82 X-l we followed 
IStrohmaver & Mushotzkvl d2003l) and extracted an 18" region cen- 
tred on the co-ordinates for the point source identified clearly in 
Chandra observations. Analysis on power spectra taken from re- 
gions s urrounding the bright ULX in the XMM-Newton observa- 
tions bv lFeng & K aaret (2007) has suggested that this is where the 
QPO originates. NGC 3628 X-l was extracted in a 23" region cen- 
tred on the co-ordinates given in Table [T] The second observation 
of NGC 55 ULX (obs. Id. 0028740101) is off-axis, causing a lower 
average count rate than observed in the first observation (obs. Id. 
0028740201) and so the two observations were treated separately. 

In order to remove background flares a light curve was ex- 
tracted from the entire chip in the 10-15 keV band and then criteria 
applied to exclude regions where the count rate stayed continuously 
over three times the mean. We only considered flares longer than 
30s. Any telemetry effects where the count rate dropped to zero for 
a period longer than 15 s were replaced with the mean count rate 
in the source file (this had only a small effect on the overall light 
curve and subsequent temporal analysis). The second half of the 
light curve from Ho IX X-l (obs Id. 0200980101) was particularly 
affected by this. The continuous regions not affected by flaring were 
defined as good time intervals. 

Luminosities were measured through fitting each of the en- 



ergy spectra with a multicoloured disc + power law model and are 
included in TableQ] 



3.1 Timing Analysis 

In order to estimate the PSD from an observation the available 
good time was split into at least five segments of equal length. 
The FFTs from each segment were averaged and the result re- 
binned in frequency giving a factor 1.1 increase in frequency per 
bin (i.e binning over / — » 1.1/), but ensuring there were at least 
20 points in each. This averaging process permits minimum-^ 2 
fitting to be used as a maximum lik elihood method dvan der Klisl 
1 989aj; iPapadakis & Lawrencd[r993h. The PSD was normalised to 
the fractional rms squared l van der Klisl 1 19971) . Visual examina- 
tion reveal ed no QPOs in addition to thos e already reported in 
M82 X-l dStrohmaver & Mushotzkvl [2003b and NGC 5408 X-l 
dStrohmaver et alj2007h . Figure[T]shows a selection of eight of the 
sample of PSDs, these are all plotted with the same segment size 
for ease of comparison, it must be noted that some of the PSDs 
used for the analysis were produced with larger segments and this 
allows a few more points of data to be visible at low frequencies. 

At high frequencies the PSD is expected to be dominated by 
Poisson noise, random fluctuations in the photon counts even in 
the absence of intrinsic variability in the source flux. T his produces 
a cons tant (i.e. flat) additional component in the PSD dvan der Klisl 
1989a) with an expected noise level (in the absence of detector non- 
linearities) of approximately 2(5 + B)/S 2 where 5 is the source 
count rate and B is the backgrou nd count rate, calculated from the 
count s in the good time intervals (Utt lev et alj200 2; Vaugh arTet al.l 
2003). In order to determine the Poisson noise level in each obser- 
vation the high frequency region of the PSD (> 0.5 Hz) was fitted 
with a constant and the normalisation compared to the expected 
Poisson noise level. Generally these correspond within the ~ 2% 
errors on the PSD estimate. However for the two brightest sources, 
Ho II X-l and M82 X-l, the measured values were ~ 3% below 
the prediction levels. These differences may be due to instrumen- 
tal eff ects such as event pile-up or dead time (see also lGoad et al.l 
120061) . 



4 RESULTS 

The following subsections describe the results of fitting various 
simple models to the power spectral data. 



4.1 Variable Sources 

In the absence of any intrinsic variability from the source the power 
spectrum is expected to be constant with frequency resulting from 
Poisson noise in the photon counting detector. A reasonable test for 
the presence of intrinsic variability is therefore a x 1 goodness-of- 
fit test comparing the data to a constant. The results are given in 
Table [2] Sources with a p-value < 0.05 were considered to show 
significant intrinsic variability (with this threshold one might ex- 
pect ~ 1 false detection in a sample of 19 observations). In fact 
six sources (eight observations) showed significant variability: Ho 
IX X-l, NGC 5408 X-l, NGC 1313 X-l, NGC 1313 X-2, NGC 
55 ULX and M82 X-l. Investigation of the PSD from NGC 4559 
X-l revealed it t o be flat with no indi cation of the variability fea- 
tures analysed in ICropper et al. (2004), this is believed to be due to 
artificial variability c aused by the analysis method as discussed in 
iBarnardet ai1d2007h . 
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Figure 1. Power spectra of 8 observations from the sample, representative of the range of variability and power seen. All light curves are split into at least 5 
equal segments in time and then re-binned in frequency to give increments of l.lfo. The Poisson noise level for each observation is also shown (dotted line). 
The top six figures show comparisons between sources which display power law type variability and sources which display no strong short term variability, 
just Poisson noise. The figure of the second observation of M82 X-1 (obs. ID. 02060801) shows variability but not of a power law type, it is visible as excess 
power in the PSD above the line indicating the Poisson noise level. The QPO in M82 X-1 at 0.12 Hz is clearly visible, as is the steep, strong red-noise feature 
at low frequencies in the PSD from the first observation of NGC 55 ULX (obs. ID. 0028740201) 
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Source 


Alternate Names 


Obs. Id. 


RA (J2000) 


Dec. (J2000) 


Duration 


Good Time 


L.v 


name 










(ks) 


(ks) 


10 39 erg/s 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


NGC 55 III X" 




002874020 1 


00- 1 5-28 9 

\JVJm 1 aJCjG . 7 


-39 13 19.1" 


34.0 


30.3 


1.21 






009X7401 01 






28 3 


26.1 


1 18 


NGC 253 PSX-2* 


NGC 253 XMMl' 


0125960101 


00:47:32.9 


-25 17 50 3' 


38.9 


31.6 


0.25 


NGC 1 31 3 X-1 c 


NGC 1 3 1 3 I IT X \J TXO 7* 


0405090101 


03- 18-198 


-66 29 09 fiJ 


121.2 


85.8 


4.25 


Mf T r 1^11 V_9 C 


NOP 1^1^ T TT YV TXO X* 


040S0Q0101 

IrrUJuyU 1U1 


(TVI 

VJ J . 1 O ..LjL.\J 


-66 36 04 3^ 


121 2 


85 8 


5 12 


NGC 2403 X- 1 d 




0164560901' 

V/ 1 VJ*T^J UU7U 1 


07-36-25 9 


+65 35 38 9' 

TUJ _U JO.? 


77.5 


58.3 


2.62 


Hnlmher? TT X-1 c 


PGC 23324 TIT XI i TXO 31* 


02004701 01 1 


08- 19-29 9 


+70 4? 1 8 6' 


79.5 


45.3 


16.6 


MR? X-1^ 


Snnrr'p "7"' 


1 1 229020 1 ' 


09-55-50 2 


+fiQ 40 47 (Y" 


27.3 


27.4 


21.0 






0206080101" 






91.0 


57.5 


20.5 


HolmberglXX-ls 


M81 X-9", PGC 28757 TJLXlA 1X0 34* 


0200980101' 


09:57:54.1 


+69 03 47.3' 


111.8 


58.7 


8.67 


NGC 3628 X-l d 


IXO 39* 


0110980101" 


11:20:15.8 


+ 13 35 13.6" 


50.3 


46.1 


4.14 


NGC 4395 X-l rf 


NGC 4395 ULX1-', IXO 53* 


0142830101 r 


12:26:01.5 


+33 31 34.7-> 


106.1 


78.6 


0.21 


NGC 4559 X-l rf 


NGC 4559 ULX2-'; ULX-7^, IXO 65* 


0152170501 1 '' 


12:35:51.6 


+27 56 01.7> 


40.2 


33.3 


9.01 


NGC 4861 ULX" 


NGC 4861 ULXlA IXO 73* 


0141150101 


12:59:02.0 


+34 51 12.0' 


15.4 


13.8 


7.90 


NGC 4945 X-2' ! 


NGC 4945 XMMl' 


0204870101 


13:05:33.3 


-49 27 36.3' 


62.5 


22.0 


0.15 


NGC 5204 X-l d 


NGC 5204 ULXlA IXO 77* 


0405690201 


13:29:38.6 


+58 25 07.5^ 


43.4 


28.5 


5.27 






0405690501 






40.0 


15.9 


5.00 


M83 ULX" 


NGC 5236 ULXy, IXO 82* 


0110910201'' 


13:37:20.2 


-29 53 47.7' 


27.0 


18.9 


0.97 


NGC 5408 X-1 1 ' 


NGC 5408 ULXIJ 


0302900101 1 ' 


14:03:19.6 


-41 22 58.7° 


130.0 


95.7 


5.28 



Table 1. Observation Details. (1) Source Name; (2) alternate names; (3) XMM — Newton observation identifier; (4)-(5) Positions from XMM — Newton 
or Chandra data; (6) XMM - Newton observation length (ks); (7) length of good time with background flares excluded (ks); (8) Observed Luminosity 
measured from energy spectral fitting (see text). Footnotes indicate references as follows: Column 1-2 footnotes indicate the identification of sources with 
the given name(s); column 3 foot notes refer to previous pow er spectral analy ses of the observations; column 5 footnotes pr ovide references to sou rce po - 
sitions, feuainazzi et alj feOOOl): ^Roberts & Warwickl feOOOl): IStobbart et alj feOOft; ^Humphrey et alj J2003h: feoria et alj J2006t): ^Colbert et all |l995h : 
^Miller et all J2004|); foewangan et alj <2004l); IStrohmaver & Mushotzkvl J2003I); iLiu & Bregmanl 12005T): ^Colbert & Ptakl 120021) ; I Winter et all koprj) ; 
iMatsumoto et alj l200lh: iFabbianol <1989l): ^Kaaret et alj fe003h: ^Cropper et alj 12004T) fcoad et alj 120061) TFeng & Kaarell J2005h iBarnard et all J2007h 
iDewangan et alj bOOeMPewangan et alj 120061) IStrohmaver et alj 120071) 



Source 



Name 


x 2 /dof 


PrniH 


Var. 


(1) 


(2) 


(3) 


(4) 


NGC 55 ULX 


128.9/63 


1.35X10- 3 


P 




155.4/63 


< xlO" 6 


P 


NGC 253 PSX2 


74.8/63 


0.255 


N 


NGC 1313 X-1 


110.0/63 


3.07xl0~ 3 


P* 


NGC 1313 X-2 


128.6/91 


1.13xl0~ 2 


P 


NGC 2403 X-1 


58.9/63 


0.193 


N 


HolmbergIIX-1 


85.3/63 


0.554 


N 


M82 X-1 


108.3/63 


3.78xl0~ 2 


Pq 




160.3/63 


< xlO" 6 


Pq* 


Holmberg IX X-1 


261.6/63 


< xlO- 6 


P 


NGC 3628 X-1 


84.7/64 


0.134 


N 


NGC 4395 X-1 


78.1/77 


0.413 


N 


NGC 4559 X-1 


78.3/63 


0.409 


N 


NGC 4861 ULX 


64.6/63 


0.375 


N 


NGC 4945 XMMl 


61.1/63 


0.746 


N 


NGC 5204 X-1 


47.2/63 


0.875 


N 




68.8/63 


0.276 


N 


M83ULX 


158.3/63 


0.068 


N 


NGC 5408 X-1 


437.9/63 


< xl0~ 6 


Pq 



Table 2. Results on fitting a constant to determine variability. (1) Source 
Name; (2)% 2 and dof values for PSD after fitting with a constant model; (3) 
Null Hypothesis Probability for PSD fit with a constant (4) Assessment of 
Variability; N denotes no observed variability, P power law type variability, 
* indicates the variability is of a white noise form, q PSD contains a QPO 



With the exception of the second observation of M82 X-1 and 
the observation of NGC 1313 X-1, which display variability in the 
powe r spectrum close to the Poisson noise level ( IDewangan et all 
l2006tlMucciarelli et alj2006h . the variable objects all show a clear 
rise towards low frequencies compared to the observations with flat 
PSDs (e.g. Ho II X-1, NGC 5204 X-1). The continuin g rise seen in 
XTE observations of M82 X -1 with a power law of -1 ( Kaar et et all 
|2006J ; iKaaret &~F ene 2007) supported by evidence of breaks seen 
in NGC 5408 X-1 suggests that the PSDs are best fit with a broken 
power law model over a broad Lorentzian. For this reason power 
law models were fitted to the variable sources, although the limited 
number of points at low frequencies made it difficult to constrain 
the model parameters. The results are given in Table[3] 

Ho IX X-1, NGC 1313 X-2, NGC 55 ULX and M82 X-1 
(obs. Id. 0112290201) were fitted best with a power law (plus a 
constant for the Poisson noise). NGC 5408 X-1, M82 X-1 (obs Id 
0206080101) and NGC 1313 X-1 wer e fitted better with a brok en 
power law (plus a constant) as used in IStrohmaver et ail d2007l) to 
model the PSD from NGC 5408 X-1. This is the first report of a 
broken power law PSD in NGC 1313 X-1 and in order to test the 
plausibility of the detection we have attempted to fit the PSD with 
a power law (+ constant) model, and compared the x' value from 
this fit to that of a broken-power law (+ constant) model which has 
two further free parameters. Analysing the difference between the 
two x 2 values both through direct comparison and the F-test sug- 
gest that the spectral break at 0.09 Hz is significant and necessary 
to the fit with p„ uU w 10~ 4 . Fig.[2]shows the standardised residuals 
after fitting these two models. 



In the cases of M82 X-1 and NGC 5408 X-1 the obvious QPOs 
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Frequency (Hz) 

Figure 2. Standardised residuals for Power law (top) and Broken Power law 
(bottom) fits to the PSD of NGC 1313 X-l. x 1 values for the fits are also 
included. 



were modelled with narrow Lorentzians (see below). The power 
law-like regions can all be fitted with slopes between -2 and -0.5. 
For those sources with previous analyses of their variability (M82 
X-l, NGC 5408 X-l and Ho IX X-l) the power law slopes, break 
frequencies and Lorentzian parameters found agree with the previ- 
ous results , with the except i on of the QPO identified in Holmberg 
IX X-l bv lDewangan et al.l (2006) which was not reproduced with 
high significance although a continuum PSD is visible. 



4.2 Establishing limits on the variability 

The other 11 observations show no strong evidence for intrinsic 
variability, despite a range of observation lengths and source fluxes, 
and therefore sensitivity to variability. In these cases th e variability 
amplit ude was constrained, following the procedure of iGoad et al.l 
(2006), by assuming the intrinsic power spectrum has one of two 
generic forms common to the PSDs of BHBs and AGN. 

The first model is a broken power law model. It is well es- 
tablished that BHBs and AGN show PSDs that are usually domi- 
nated by a broad-band noise component with a form that resembles 
to first-order a broken power law. Most sources show a frequency 
range in which the spectral index (PSD slope) is appro ximately — 1 
- noise with this PSD are known as "flicker noise' ' jPresslll978l) . 
The integral of a spectrum diverges (if rather slowly) to low 
and high frequencies, and so must break to flatter and steeper slopes 
respectively. Therefore, we may use a broken power law as an ap- 
proximate model of the generic broad-band noise PSD. 

In many AGN and "high/soft" state BHBs the low frequency 
break (to an index flatter than -1) is not observed and must occur 
at very low frequencies iReig et al. 12003; Mark owitz et alj|2003l ; 
lUttlev et alj 12002: IReig et alj 112003: IChurazov et alj 1200 lh . This 
means a good approximation to the observed PSDs can be made 
by considering a singly broken power law (BPL) with indicies of 
-1 and -2 below and above some break frequency, respectively. 



P(f) : 




+ Cp where f < fb 
+ Cp where f > fb 



(1) 



where A is a normalisation, fb is the (high frequency) break fre- 
quency, and Cp is a constant to account for the Poisson noise. 

If we represent this spectrum in / x P(f) form the / _1 part 



of the intrinsic (i.e. after subtracting the contribution due to Pois- 
son noise) spectrum (/ < fb) becomes flat, with a constant level 
Cyf - fbP(fb) = A/t>- ^ h as been suggested that the value of 
Ci/f is fairly consta nt across BHBs and AGN, with v alues in the 
range 0.005 - 0.03 dGoad et al.ll2TO6l : |Papadakisll2004l) . Of course 
the break frequencies /j, vary widely from source to source - from 
10~ 7 Hz for AGN with BH masses ~ 10 8 M Q , to ~ 10 Hz in 
stellar mass BHBs - and on long timescales the break frequency 
may vary by as much as a decade in frequency for a single sourc e 
l lBelloni & Hasinger|[l990l ; IUttlev et al1l2002l : lBelloni et alj|2002l) . 
For the non- varying ULX observations there was generally no ob- 
servational constraint on f, and so C^f was estimated (or con- 
strained) as a function of break frequency over the range 10~ 4 - 1 
Hz. Figure[3](left panel) shows a contour plot of the upper limit on 
C Uf as a function off, for NGC 5204 X-l (obs. Id. 0405690201). 
The contours show the value of Ci/f at which the^- 2 fit statistic in- 
creased by 1.0, 2.71, and 6.63 over the minimum. These correspond 
to the 68.3, 90 and 99 per cent quantiles of the^ 2 distribution with 
y = 1 degrees of freedom, and are conventionally used to delimit 
confidence regions on one parameter. Here they are used to define 
approximate limits on the magnitude of Ci/f as a function of f b . Ta- 
MeSgives 90% limits on the value of Ci/f for each of the sources in 
the sample calculated assuming two different values of /(,, namely 
10" 3 Hz and 1 Hz. 

For the variable sources for which an unbroken power law was 
a good fit we estimate the value of Ci// by fitting a BPL model. 
Where the observed power law index was w -2 (i.e. in NGC 55 
ULX and NGC 1313 X-2) it was assumed that the high frequency 
break, fb, must lie below the observed frequency range (< few x 
10~ 3 Hz), and so we assumed fb = 10 -4 Hz. If, on the other hand, 
the observed power law index was close to -1 (i.e. Holmberg IX 
X-l, and the first observation of M82 X-l), it was assumed that the 
fb must lie close to, or higher than, the highest frequency at which 
ULX variability can be detected above the Poisson noise. In these 
cases we assumed /b = 1 Hz. 

In the cases of NGC 1313 X-l, NGC 5408 X-l and the sec- 
ond observation of M82 X-l, where a break in the power law is 
detected, the observed indices at low frequencies are rather low 
(closer to than -1) indicating that the observed PSD covers the 
low-intermediate part of the model. Here it is flat, breaking from 
an index of to -1 the low frequency flattening discussed earlier. 
In these cases the intermediate-high part of the PSD was again as- 
sumed to break from an index of -1 to -2 at f b = 1 Hz, while 
the low frequency index of was fixed, the location of the lower 
frequency break was left as a free parameter. 

The second model represents the band-limite d noise (BLN) 
commonly observ ed in BHB low/hard states dNowakl [2000; 
lBellonietai]|2002h . In this model the variab ility is mostly lim - 
ited to ~ 1-2 decades of temporal frequency ( Ivan der Kli j|2006l) . 
A simple model for this type of noise is a broad (incoherent) 
Lorentzian: 



P(f) 



R 2 Q.fo/n 



/ 2 + G 2 (/-/o) 2 



+ C P 



(2) 



where R is a normalisation term (approximately integrated power, 
lNowakll2000l) . fo is the resonant frequency, and Q is the qual- 
ity factor. Again, C P is a constant to account for the Poisson 
noise. The Q-factor defines the frequency width of the Lorentzian, 
Q x f /HWHM where /o is the central Lorentzian frequency and 
HWHM its half width at half maximum. Typically a feature with 
i ? > 2 is classified as a Q PO and Q < 2 as band-limited noise 
i van der Kli j[l989al. [2006). In order to model band-limited noise 
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Source a[ f\, (10 2 Hz) a<i A\ C p x 1 /dof 

(1) (2) (3) (4) (5) (6) (7) 



NGC 55 ULX (1) 


1.75±0.2 






(1.53±0.6)xl0 -5 


1 7 1 +0.06 
-0.02 


86.9/70 


NGC 55 ULX (2) 


1.96±0.04 






(6.55±2.0)xl0" 6 


3.87±0.02 


122.6/70 


NGC 1313 X-l 


-0.082±0.07 


Q l +0.6 
*- l -0.9 


2.35±1.0 


0.26±0.06 


3.16±0.008 


71.9/70 


NGC 1313 X-2 


1 Q 1+0.05 
L - y -0.1 






(5.82±3.0)xl0~ 7 


3.59±0.007 


69.3/70 


M82 X-l (1) 


().70±0.()6 






(l.l±0.3)x 10~ 2 


1.346±0.006 


41.6/55 


M82X-1 (2) 


14 +a02 

-0.05 


2.0±1 




(5.5±2.2) x 10~ 2 


0.847±0.002 


66.7/70 


HolmbergIXX-1 


0.48±0.04 






(2.55±0.3)xl0~ 2 


1.36±0.04 


73.6/70 


NGC 5408 X-l 


-0 10+ 007 
u.iu _ 003 


n n+OM 
"-0.07 




4.0±0.8 


2.31+0.01 


100.2/77 



Table 3. Best fit models for the sample sources displaying variability. (1) Source names; (2) Power law index (below break for broken power law); (3) Break 
frequency of broken power law component in Hz; (4) power law index above the break; (5) Normalisation of power law/broken power law component; (6) 
Normalisation of Poisson noise component; (7)^ 2 and dof values for the fitted models. A L orentzian component (with Q = 5) w as included in the model s 
fitted to the QPOs in M82 X- 1 and NGC 5408. For more detailed examination of the QPOs see lStrohmaver & Mushotzky 1 2003) and Stroh maver et all fc007l) . 




^O- 3 0.01 0.1 1 : : 

Break Frequency f b (Hz) 10 -t 10 -> 0.01 0.1 1 10 10-' 10~ 3 0.01 0.1 1 10 

Central Lorentzian Frequency (Hz) Central Lorentzian Frequency (Hz) 



Figure 3. Contour plots for our three variability models from NGC 5204 (obs Id. 0405690201), limits shown correspond to Ax = 1.0, 2.71 and 6.63. Left: 
High/soft model, broken power law -1/-2. Centre: Band-limited noise, Lorentzian (Q=0.5). Right: Narrow Lorentzian (Q=5.0) 



a value of Q = 0.5 was assumed. The same model form but with 
Q = 5.0 was used to model QPO features. Of course, given the 
quality of the data the resonant frequency /o of the BLN could not 
be constrained, and so the strength of any BLN was estimated (or 
constrained) as a function of f spanning 10~ 4 - 1 Hz, in the same 
manner as with the BPL model. In those observations that showed 
no obvious QPOs the same procedure was used to put limits on the 
strength of any QPOs as a function of QPO frequency. The centre 
and right panels of Fig. [3] show the upper limits on amplitudes of 
the BLN and QPO models for NGC 5204 X-l (obs Id. 040569201). 
These were found through minimising the fit to the models over 
1000 different values of either the central Lorentzian frequency or 
the break frequency over the XMM-Newton frequency range. 

Table [4] gives the approximate limits on the amplitudes of the 
two models (BLN and BPL) assuming two representative values of 
the characteristic frequency in each case (1 mHz and 1 Hz). The 
limits were calculated using a A^ 2 = 2.71 criterion (which corre- 
sponds to the 90 per cent limit of a^ 2 distribution with v = 1). The 
fainter sources (such as NGC 4945 X-2, NGC 4395 X-l, NGC 4861 
ULX, NGC 3628 X-l and M83 ULX) have upper limits such that 
variability from these sources would have to be particularly strong 
to have been detectable. For example, in NGC 4395 X-l a BLN 
spectrum could have R x 20% fractional rms and this would barely 
reach the 90% confidence limit even assuming a characteristic fre- 
quency best-suited to detection (f ~ 10~ 3 Hz). For those sources 
in which variability was detected the measured values of the am- 
plitude are given for each model. The observation of NGC 5408 



X-l gives the amplitude of a BLN spectrum (assuming f & 10 3 
Hz) as R x 18% whereas the upper limit on the amplitude of the 
same model in the case of NGC 5204 X-l (obs Id. 0405690201) 
gives R < 7% even though the two observations show comparable 
Poisson noise levels. 



4.3 Limits on QPOs 

The data were examined for possible QPOs by re-binning the FFT 
data in both frequency and time using various different bin widths. 
As discussed above a narrow (Q = 5) Lorentzian was used to model 
a QPO. For each observation the change in^- 2 was calculated as a 
function of the two free parameters R 2 and f . An improvement 
of A^ 2 = 11.3 was considered indicative of a QPO. This process 
clearly recovered the known QPOs in NGC 5408 X-l and M82 X- 
1 but did not find any further candidates. The known QPOs show 
typically R 2 ~ 10% and f ~ 2 - 5 x 10~ 2 Hz. By comparison the 
limits on the QPO strength derived from observations of several 
other bright ULXs, such as NGC 5204 X-l and Ho II X-l, show 
that such QPOs would have been detected i f they were present (see 
Fig.[D. The QPO claimed in Ho IX X-l bv lDewangan et alj < l2006|) 
was not reproduced with high significance. In these sources QPOs 
similar to those in NGC 5408 X-l and M82 X-l sources would 
have to be fairly muted if present. 
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Source 


Broken Power Law (BPL) 


Band limited noise (BLN) 


Name 


Ci// (icr 4 ) 


C Uf (in- 4 ) 


R 2 (10" 2 ) 


R 2 (10~ 2 ) 


R 2 (10- 2 ) 




10^ Hz. 


IHz, 


\Q- 3 Hz 


I Hz 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


NGC 55 TIT X 


- 
- 


177±46 
463±142 


- 


• 


3.2±0.5 
9.2±2 


NGC 253 PSX2 


< 65.4 


< 18.5 


< 1.64 


< 35.2 


• 


Nnr 1 3 1 3 x- 1 


- 


136±46 


- 


- 


6.0+? 8 


Nnr mi x-2 

1 1 V J v . T J T J yv 




26±14 


- 




°- 65 -o:l 


Nnr 24oi X- 1 


< 2.9 


< 2.3 


< 0.282 


< 15.5 


- 


Hnlmherp TT X-1 


< 2.88 


< 3.7 


< 0.0610 


< 2.23 




M82 X-1 




7.4±5.0 
22.1 ±6.0 






4.2±2 
2.5±0.5 


HolmbergIXX-1 




6.0±1.7 






0.88±0.2 


NGC 3628 X-1 


< 170 


< 43.1 


< 4.49 


< 33.1 




NGC 4395 X-1 


< 260 


< 64.7 


< 5.58 


< 92.8 




NGC 4559 X-1 


< 7.11 


< 28.4 


< 1.72 


<43.4 




NGC 4861 ULX 


< 401 


< 99 


< 10.8 


< 206 




NGC 4945 XMM1 


< 180 


< 51.8 


< 5.37 


< 293 




NGC 5204 X-1 


< 21 


< 5.9 


< 0.494 


< 10.5 






< 34.3 


< 5.4 


< 0.539 


< 13.3 




M83 ULX 


< 29.1 


< 64 


< 7.56 


< 42.3 




NGC 5408 X-1 




104±11 






2.99±0.5 



Table 4. Upper limits on the fitted model amplit udes using t wo different power spectral models. The estimate of Ci/y from the BPL model at f;, = 10 3 Hz for 
Ho II X-1 differs from the value found in lGoad et al " (2006) differs due to the use of different binning for the fitting process. Columns: (1) Source Name; (2) 
Ci// from the broken power law model with f/, = 10~ 3 Hz; (3) For non-varying sources: Ci// from the broken power model with the fj, = \Hz\ For varying 
sources: Ci// from fit to model according to shape of PSD: See text; (4) Upper limit on fractional rms within a Lorentzian of coherence Q=0.5 fit at 10~ 3 Hz; 
(5) Upper limit on fractional rms within a Lorentzian of coherence Q=0.5 fit at 1 Hz; (6) Fit normalisation values to Band limited noise when present. All 
limits correspond to Ay 2 = 2.71 (which is approximately 90% confidence). 



5 DISCUSSION 

The main results of the temporal analysis of 19 observations of 16 
sources may be summarised as follows: 

• Six sources (eight observations) showed significant intrinsic 
flux variability in the lfr 4 - 1 Hz range (NGC 5408 X-1, NGC 
1313 X-1, NGC 1313 X-2, NGC 55 ULX, M82 X-1 and Ho IX 
X-1). 

• Of these six, two show strong QPO features (NGC 5408 X-1 
and M82 X-1, as previously identified) and all show a continuum 
spectrum consistent with a power law or broken power law model. 

• Of the remaining 10 sources, the intrinsic variability ampli- 
tudes were constrained by comparing the data to simple PSD mod- 
els based on BHB observations. 

• The limits obtained for several of the observations (e.g. NGC 
4559 X-1, NGC 5204 X-1, Ho II X-1) are such that the strength 
of the intrinsic variability, in the frequency range observed, must 
be substantially lower than the observed power seen in the variable 
sources. 



5.1 Break Frequencies and BH mass 

It is in principle possible to draw inferences about the nature of 
the putative accreting black hole in the ULXs if analogies can 
be formed between their observed properties and those of BHBs 
and AGN. Perhaps the simplest and best studied relation linking 
BHBs and AGN is the correlation between the PSD break frequen- 
cies and BH mass (e.g. Markowit z~et al.ld2003 l): lDone & Gierlinski 
(2005); McHardv et a l. (2006)), a correlation which may well in- 
clude ULXs. 



The 3 PSDs of NGC 5408 X-1, M82 X-1 (obs. ID. 
0206080101) and NGC 1313 X-1 show a break from an index 
of ~ -1 above the break to ~ below, which resembles the 
Tow frequency break' observed in low/hard BHB observations 
dBelloni & Hasingeilll990l : lNowakll200d iMcClintock & Remillarj 
l200q) . ISoria et al.l ( 120041) have used the break frequency in NGC 
5 408 X-1 to derive an exp ected BH mass of ~ 100 M e . For M82 X- 
l lDewangan et al.l d2006l) measured the PSD break frequency in re- 
lation to the QPO frequency and drew a comparison with BHB fea- 
tures i n order to estimate the BH mass at 25-520 M Q . Casell aet al.l 
d2008l) applied the relatio nship in Kordi ng et al] ( 12007^ which ex- 
tends the analysis from IfvIcHardv et al.l d2006D into systems in a 
low/hard state, to estimate a BH mass between 1 15 - 1305 M Q for 
NGC 5408 X-1 and 95 - 1260 M Q for M82 X-1. 



By contrast the break in NGC 1313 X-1 at 0.09 Hz is compa- 
rable to those observed in BHBs, for example the observed break of 
Cyg X-1 in the low/hard state is 0.08 Hz. Cyg X-1 is thought to have 
a BH mass of ~ 10M o . However the observed luminosity of Cyg 
X-1 in this state is only ~ 4x 10 37 erg/s dNowak et al.ll999h a factor 
of 100 lower. Either NGC 1313 X-1 hosts a commensurately larger 
BH mass than Cyg X-1, or it does not but is accreting at a much 
higher relative rate (L/L Edd > 1), or at least appears to be due to 
anisotropic emission (or some combination of these). In either case 
the flat PSD below 0.09 Hz is difficult to account for. The flat slope 
is most commonly observed in Tow/hard' state BHBs, but this state 
is tho ught to occur when L/Lem < 0.05 dMcClintock & Remillard 
120061) . Assuming NGC 1313 X-1 to be in such a state would there- 
fore require M BH > 600M o , yet the frequency of the break matches 
well the expectation for a Mbh ~ 10M Q black hole. This supports 
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the s uggestion that U LXs may represent a different, 'ultraluminous' 
state (Roberts 2007). 

NGC 1313 X-2, NGC 55 ULX and Ho IX X-l show excess 
power rising to low frequencies below ~ 5 x 10~ 3 Hz. As can be 
seen in Figure [T] the PSD from the first observation of NGC 55 
ULX (obs. ID. 0028740201) shows a steep red noise region at low 
frequencies, also visible in the PSD from the second observation 
(obs. ID. 0028740101). When the two observations are fitted with 
a power law + constant model, spectral indices of 1.75±0.2 and 
1.96±0.4 are found respectively, with no sign of a low frequency 
break. This would suggest, assuming the system is in the high/soft 
state, that the high frequency break would have to be bel ow 10~ 4 
Hz, in dicative of a BH with a mass in excess of 10 6 M Q dRobertsI 
l2007t> more typical of AGN. However it mu st be pointed out tha t 
dips in the light-curve have been identified dStobbart et al]l 2004) 
and that these features affect the PSD on timescales commensurate 
with the strong variability power seen in its PSD. The power laws 
seen in the PSDs from Ho IX X-l and NGC 1313 X-2 have spectral 
indices of ~0.5 and ~ 1 .9 1 respectively. The few rising points visi- 
ble in the PSD of NGC 1313 X-2 do not allow for accurate fitting 
and so little can be derived from this value. For Ho IX X-l however, 
the slope is much clearer and we have carried out a basic analysis. 
We have been able to constrain the likelihood of a brake to a slope 
of -2 to be abov e 0.1 Hz. This would indicate a BH with a mass 
below 1000 M dRobertsll 20071) . 

The possibility of estimating the masses of ULXs by com- 
paring characteristic (e.g. break or QPO) frequencies with those 
observed in BHBs and AGN, is dependent on there being a suit- 
ably strong analogy mapping PSD characteristics between the 
source, together with a simple mass-frequency scaling relation 
(McH ardv et~al] |2006). The r ecent discovery of an apparent high 
frequency QPO within a AGN dGierliriski et al. 2008) indicates that 
these features are common across the range of observed sources. If 
however ULXs do not display variability in forms typical to either 
BHBs or AGN and instead exist in a new 'Ultraluminous state' 
dRobertsI [2007h with distinct and different characteristics to those 
states observed in AGN and BHBs, then these mass estimates will 
no longer be reliable. 

5.2 Variability vs Flux and Luminosity 

The two main properties of an observation that determine the sensi- 
tivity to intrinsic spectral power at a given frequency are the length 
of the observation T (the number of 'cycles' of that frequency 
that were observed) and the count rate of the source (which de- 
termines the Poisson noise in the observation). Denoting the mean 
count rate as / we may approximate the significance (in units of 
<x) of an excess in power over the Poisson no ise level by n^if) « 
i IPs (/) yr x A/(/) (see Ivan der Klisl["l989bl) where Ps(f) is the 
power density due to the intrinsic variations in the source flux, and 
A/(/) gives the frequency binning. Thus, with similar observation 
lengths T and frequency re-binning A/(/) the sensitivity should 
scale as n& oc IPs(f) an d so if a detection is made when exceeds 
some number the limit on the power that can be detected scales as 

p s (f)^r l . 

Fig. E] (left panel) shows the estimated variability amplitudes 
(detections and upper limits) against source flux. Also on these 
plots are typical values for Ci/; observed in a BHB Cyg X-l and the 
AGN Mrk 335, illustrating the consistency in Ci/f between BHBs 
and AGN. The upper limits clearly improve (i.e. give tighter limits) 
at higher fluxes, as expected. In particular, they indicate it is in gen- 
eral not possible to detect even strong variability (Ci/f > 0.05) in 



sources fainter than ~ 10~ 12 erg s cnT 1 . Fig.|4](right panel) shows 
the same variability amplitudes against observed X-ray luminosity, 
with no strong trend apparent. Above 10 3 ' erg s' 1 (i.e. definite ULX 
luminosities) there is a range of ~ 30 in the variability amplitude, 
confounding any efforts to identify source states on the basis of 
luminosity-variability correlations. Considering only those source 
for which variability was detected (filled circles) there appears to 
be a correlation suggesting a decrease in variability with increas- 
ing luminosity. However there are two reasons why this could be 
misleading. The first is that the two highest estimated values of 
Ci/f are the observations of NG C 55 ULX. This sour ce shows dip- 
ping episodes in its variability dStobbart et al] |2004). not seen in 
other ULXs, and this anomalous behaviour enhances the variabil- 
ity amplitude. Removing these points substantially decreases the 
apparent significance of any luminosity-amplitude anti-correlation. 
Secondly, including the limits from observations that did not have 
variability detections shows that some lower luminosity sources 
(~ 10 39 erg s" 1 ) must have very low amplitude variability, opposite 
to the prediction of a general luminosity-amplitude anti-correlation. 
When these points are considered there appears to be no obvious 
correlations within the plot. 



5.3 Are some ULXs significantly less variable? 

The upper limit contours from four bright sources (Ho II X- 1 , NGC 
2403 X-l, NGC 5204 X-l and NGC 4559 X-l) show that these 
sources are significantly less variable than others (such as NGC 55 
ULX, NGC 5408 X-l, etc.). Fig. [5] compares the data from NGC 
5204 X-l with the model obtained by fitting the data from NGC 
5408 X- 1 (these two have similar luminosities and count rates, and 
hence Poisson noise levels in their PSDs). The figure clearly shows 
the variability seen in NGC 5408 X- 1 to be absent or greatly sup- 
pressed in NGC 5204 X-l. 

The difference between the strength of variability observed 
in many ULXs and that expected by analogy with BHBs was 
confirmed using simulated data. Two PSDs of Cyg X-l - one 
showing BLN ("low/hard state") and one showing BPL ("high/soft 
state"Jj - were used as examples of BHB variability. Random 
time series were genera ted from these PSDs (using the method of 
iTimmer & K onig 1995) and given the same mean count rate and 
exposure time as the XMM-Newton observation of Holmberg II X- 
1. Poisson noise was included at the appropriate level, and the sim- 
ulated data were analysed in the same manner as the real ULX data. 
Variability was clearly detected (i.e. a constant was rejected by a^ 2 
goodness-of-fit test) in the simulated data for both PSD shapes. The 
same was true when the PSD frequencies were reduced (by factors 
of up to 10 5 ) prior to generating the time series - i.e. simulating 
the effects of larger black hole masses by scaling the variability 
timescales. 

The results confirm that if Holmberg II X-l was variable with 
a PSD shape and amplitude similar to that seen in Cyg X-l its 
variability should have been detected (at high significance) in the 
XMM-Newton observations. This is true whether the PSD shape is 
BPL or BLN, and holds when the characteristic variability frequen- 
cies are reduced as might be expected for IMBHs. The results of the 
simulations confirm the the PSD fitting results (discussed above), 
namely that the apparent lack of variability in the observations of 



3 The observations were taken from the RXTE archive: the BLN data came 
from proposal P10236 made on 1996 December 16-18 and the BPL data 
came from P10512 made on taken 1996 June 18. 
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Figure 4. Upper limits to high/soft model against Flux (left) and luminosity (right) for fi, = 1 mHz. Points de fined as va rying are represented b y filled circles, 
whils t non-varying sources are unfilled. Also shown are the measured values of Ci// for BHB Cygnus X-l l Uttlev 2007) and the AGN Mrk 335 iArevalo et al. 
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Figure 5. Average power spectra of NGC 5204 (obs Id. 0405690501) with 
the best fitting model for NGC 5408 X- 1 , demonstrating clear absence of 
variability in the former source even though the observed luminosities from 
the two observations are very similar. 



sources such as Holmberg II X-l and NGC 5204 X-l must be in- 
trinsic to the sources, not a result of insufficient observations. 



5.4 Interpreting the lack of variability 

There are two fundamentally different ways to account for the lack 
of variability in the ULX observations: as an observational effect 
caused by the limited energy or frequency bandpass of the observa- 
tions, or as a result of the sources being intrinsically under- variable. 
We consider some specific possibilities below. 

The energy spectra of BHBs in their higher accretion rate 
states are often dominated by a (quasi-)thermal component as- 
sociated with the accretion disc. By comparison with the non- 
thermal (power law-like) hard X-ray contin uum, this spectral com - 
ponent shows little short term variability dChurazov etai]|200lh . 
LMC X-3 is an example of a BHB that often shows very low am- 
plitude short term variabil ity in the "high/so ft" state (dominated 
by thermal disc emission; iNowak etalj|200ll) . See also a discus- 



sion of this point in the context of BHB-AGN comparisons in 
iDone & Gierlinskil §005). It could be that a similar component 
dominates the spectra of the under- variable ULXs within the XMM- 
Newton bandpass. If this was the case we would expect the thermal 
emission to peak around ~ 1 — 2 keV, higher than expected for 
IMBH models, which predict a much cooler te mperature ~ 0.2 
keV dMiller et alj[2003l |2004 IStobbart et alj|2006l) , i.e. this solu- 
tion requires BH masses similar to Galactic BHBs. One problem 
with this hypothesis is that the energy spectra of ULXs often appear 
more complex than the thermal-dominated spectra of "high/soft" 
BHBs when viewed in the hi g hest quality XMM-New ton obser- 
vations I Stobba rt et al. | d2006h : iGladstone et al. (2009)), although 
even then there are som e ULXs that do appear t o have modi- 
fied di sc-like spectra (e.g. Ivierdavanti 'et alj l l200r3) : iHui &~K rolik 
(2008)). A prediction of this hypothesis is that the sub-variable 
ULXs should show stronger variability at higher energies, where 
the sp ectrum is domina ted by variable non-thermal emission. 

iGoad et al.1 (2006) suggested that the lack of variability seen in 
Holmberg II X-l may be related to the so-called '^-state' observed 
in GRS 1915+105, where the variability power is concentrated in 
a BLN spectrum at high frequencies (> 0.1 Hz) and the re is often 
little variability on longer timescales (B elloni et alj200oT) . The inte- 
grated power over the 10 -4 - 0. 1 Hz range is typica lly only ~ 10~ 3 
(fractional variance; see e.g. Zdzi arskiet"ai] [2005) which is close 
to the detection limits of the ULX observations. Again, this expla- 
nation requ ires BH masses com parable to that of GRS 1915+105 
(~ 14 M ; Greineret al. 2001) as a significantly larger BH mass 
would shift the BLN variability into the frequency bandpass of the 
ULX observations (10~ 4 - 0.1 Hz) where it would be detectable. A 
prediction of this hypothesis is that ULXs should show significant 
variability power at higher frequencies (e.g. 0.2 - 20 Hz), in the 
form of a BLN component and possibly strong QPOs, as does GRS 
1915+105. 

The above explanations hypothesise that the variability in 
ULXs is similar to that of BHBs (and AGN) but is often not ob- 
served due to the limited energy and/or frequency bandpass of the 
available data. Another possibility is that the X-ray emission from 
ULXs is intrinsically less variable (on frequencies > 10~ 3 Hz) than 
expected by analogy with other accreting systems. 

The timescales involved are sufficiently large to rule out at 
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least one model for the suppression of variability within the ULX 
system: scattering in an ext ended region. Recent studies of the 
X-ray spectra of ULXs (e.g. bone & Kubotall2006l : IStobbart et al.1 
l200alGoadeTai]|2006h have put forward the idea of an optically 
thick corona of hot electrons that modifies the emerging energy 
spectrum. The timescale for the diffusion of photons through an 
electron scattering medium is of order At ~ Rr es /c, where R is 
the size of the scattering regi on and T es (> 1) is t he optical depth 
to ele ctron scattering (see e.g. lNowak et alJ["l999l : lMivamoto et al.l 
1 1988L and references therein). On timescales shorter than At the 
variability will be strongly attenuated. The ULX observations show 
weak or absent variability on timescales at least as large as ~ 10 2 
s, which would require Rr es ~ 100 lt-s. This is implausibly large; 
even with Mbh ~ 100 M G an optical depth of t ei ~ 10 - 20 (see 
Gladst one et alj|2009h would need a scattering region as large as 
R ~ 2 x 10 9 m (~ 7 x 10 3 gravitational radii), comparable to the bi- 
nary separation of Cyg X-l. This analysis suggests electron scatter- 
ing is not a viable mechanism for suppressing variability in ULXs 
on such long timescales, although this is still a promising model 
for the energy spectra which require a much smaller emission re- 
gion than needed here. 

Of course, it remains plausible that ULXs are examples of 
an accretion state that is simply not found (or very rare) among 
the commonly observed BHBs and AGN, and which displays 
very stable X-ray emission. Nevertheless, any mechanism for ex- 
plaining the lack of variability in some ULXs faces a challenge: 
the X-ray spectra of the variable (e.g. NGC 54 08 X-l) and non- 
variable (NGC 5204 X-l) ULXs are quite s imilar JSoria et alJiooH ; 
iRoberts et all 12004 IStobbart et ail [2006) while their variability 
properties appear very different (see Fig. [5}. This apparent incon- 
sistency requires further study to understand fully, multiple obser- 
vations of a few ULXs could allow for further clarification of this 
problem. 

Longer observations with XMM-Newton would allow the fre- 
quency range to be extended to lower frequencies and could re- 
veal in more detail the PSDs of ULXs such as NGC 55 ULX and 
NGC 1313 X-2, allowing better comparisons with BHBs and AGN. 
Observations of ULXs with more powerful X-ray missions (such 
as Astro-H and the proposed International X-ray Observatory) are 
needed before ULX variability can be studied at the higher frequen- 
cies and in the higher energy bands needed to test the hypotheses 
discussed above. 
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